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Abstract. We present and study in this paper a simple algorithm that produces 
so called growing Parallel Random Apollonian Networks (P-RAN) in any dimension 
d. Analytical derivations show that these networks still exhibit small-word and scale- 
free characteristics. To characterize further the structure of P-RAN, we introduce 
new parameters that we refer to as the parallel degree and the parallel coefficient, 
that determine locally and in average the number of vertices inside the (d+l)-cliques 
composing the network. We provide analytical derivations for the computation of the 
degree and parallel degree distributions, parallel and clustering coefficients. We give 
an upper bound for the average path lengths for P-RAN and finally show that our 
derivations are in very good agreement with our simulations. 
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1. Introduction 

During the last decade, the study of network topologies has become a useful way to tackle 
the understanding of information flow within complex natural or artificial systems. The 
applications range from sociology, logistics, epidemiology, immunology, neural networks 
characterization, granular packing analysis, networking, etc. Among a multitude of 
proposed models, scale-free and small world networks have been widely addressed, 
essentially because many empirical or real life networks display such properties P, [9]. 
This is the case for random graphs, social networks, the web and for gene networks 
for instance. Basically, scale free networks display a power-law degree distribution, 
p(k) ~ /c~ 7 , where k is connectivity (degree) and 7 the degree exponent [2], while in 
small world networks, most vertices can be reached from any other by a small number 
of hops or steps. Small world networks are characterized by a high clustering coefficient, 
e.g. a high level of vertices interconnection, and small average path length, namely small 
minimum path length in average between any pairs of vertices in the network. 




Figure 1.1. 2D Apollonian gasket and corresponding network. 1 st generation: disks, 
2 nd generation: squares, 3 rd generation: triangles 

Among the topologies that display scale free and small world properties, Apollonian 
networks [1] have recently attracted much attention [8j Ej. Apollonian networks are 
constructed from a fractal generated from a set of hyper-spheres, where any hyper- 
sphere is tangent to the others. This fractal is also known as the Apollonian gasket, 
named after Greek mathematician Apollonius of Perga. The 2D Apollonian network, or 
Deterministic Apollonian Network (DAN) [ilj, is obtained by connecting the centers of 
touching spheres (interpreted as vertices) in a three-dimensional Apollonian gasket by 
edges, as shown in Fig 11.11 The first generation for this fractal network is characterized 
by disks vertices, the second generation is characterized by square vertices and the third 
generation is characterized by triangle vertices. Extension to higher dimension have 
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been provided in [TT] . 

Ramdom Apollonian Networks (RAN) [15] , differ from the recursive construction 
of DANs, as a RAN starts from a (d+l)-clique (a triangle in dimension 2) containing 
d + 1 vertices. Then, at each time step, a single (d+l)-clique is randomly selected from 
the set of (d+l)-clique in the network that do not already contain a vertex connected 
to all the vertices composing the (d+l)-clique. The selected (d+l)-clique is then used 
to insert a new vertex linking to all of the d + 1 vertices of the selected (d+l)-clique. 
2D Random Apollonian Networks (RAN) have been extensively studied in [T5l IT3] . and 
extension to high dimension RAN (HDRAN) provided in [T2] . 

Some recent attempts to make use of RAN like structures in P2P application [refs] 
faces the requirement to maintain such topologies in dynamic conditions, e.g. when 
vertices almost freely enter and leave the network. For RAN or HDRAN topologies, 
the repairing process when vertices leave the network is quite costly and limits the 
range of potential applications. In order to simplify the topology repairing process 
(that is beyond the scope of this paper), we are considering an extension of the RAN 
or HDRAN topologies to what we call Parallel RAN (P-RAN). This new topology 
that differs slightly from RAN or HDRAN allows to insert several vertices inside a 
(d+l)-clique, each inserted vertices being fully connected to all the vertices composing 
the clique. This extension constructs parallel random Apollonian structures that we 
formally study through out the paper. 

After a short presentation of Parallel Deterministic Apollonian Networks (P-DAN) 
and Parallel Random Apollonian Networks (P-RAN) in the first two sections, we 
introduce in the third section the parallel degree distribution and parallel coefficient 
for such networks and study their asymptotic statistical properties for any dimension. 
The fourth, fifth and sixth sections give the derivations respectively for the degree 
distribution and the degree exponent, the clustering coefficient and the average path 
length for P-RAN. Extensive simulation results are provided through out these sections 
to validate as far as possible the analytical derivations. A short conclusion ends the 
paper. 

2. Parallel Deterministic Apollonian Networks 

A parallel deterministic Apollonian network in dimension d is constructed recursively 
from an initial (d+l)-clique allowing to insert at step t more than one vertex into 
(d+l)-cliques composing the network at step t — 1. Various rules can be adopted for 
the construction of Parallel Apollonian networks. Some of them lead to Expanded 
Apollonian networks [10] or recursive clique trees [4] for which at each time step, a 
new vertex is inserted in every (d+l)-clique composing the network. In the following 
subsection, as an example, we propose other rules that lead to a different topology. To 
characterize the parallel nature of this kind of networks, we introduce what we call the 
parallel degree m > of a (d+l)-cliques that characterizes the number of vertices inside 
the clique and fully connected to the vertices composing the clique. This constructing 



Parallel Random Apollonian Networks 
process is detailed in Algorithm [T] 



4 



Data: d: dimension of the P-DAN; tMax: maximum number of steps; m the 

parallel degree 
Result: r a d-dimensional P-DAN 
t *- 0; 

Initialize r to a (d+l)-clique, c (r contains 1 (d+l)-cliques); 
C <— c the set of (d+l)-cliques composing the P-DAN ; 
while t < tMax do 

for all c in C do 

Insert m new vertices into c, fully connected to the vertices composing c ; 
Insert into C the m.id + 1) new created (d+l)-cliques ; 
end 

t <-t+ 1 ; 
C^C"; 
end 

Algorithm 1: P-DAN constructing algorithm 




Figure 2.1. 2-dimensional P-DAN at t=0 (left), t=l (middle) and t=2 (right) 



Constructing algorithm 

Following Algorithm [1] specification, initially, a network containing d + 1 vertices and a 
single (d+l)-clique is created. At each time step t, m > 1 vertices are added into all 
existing (d+l)-cliques created at time step t — 1 in the current network and each new 
vertex is connected to each vertices of the embedding (d+l)-clique, creating m.(d + 1) 
new (d+l)-cliques. Figure [27T1 presents the first three steps of the P-DAN constructing 
algorithm. 

3. Parallel Random Apollonian Networks 



We define Parallel Random Apollonian Networks as RAN for which a new vertex can 
be inserted at time step t in any (d+l)-clique composing the network, whatever its 
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creation time step is. This means that a (d+l)-clique can contain in its inside more 
than one vertex fully connected to the vertices composing the clique as detailed in 
Algorithm [2j To our knowledge, no previous work have been reported specifically on 
P-RAN. Nevertheless, some similarity can be found for simple topologies described in 
one dimension in [5]. We study in the following sections P-RAN for any dimensions. 




Figure 3.1. 2-dimensionnal P-RAN. One vertex is added to a randomly chosen 3- 
clique at each time steps. Edges added at each time step are dashed 

3.1. Constructing algorithm 

Data: d: dimension of the P-RAN; tMax: maximum number of steps 
Result: r a d-dimensional P-RAN 
t «- 0; 

Initialize r to a (d+2)-clique (r contains d+ 1 (d+l)-cliques); 

C <— cl,c2,c3 the set of (d+l)-cliques composing the initial (d+2)-clique ; 

while t < tMax do 

Select randomly a (d+l)-clique, c, in C ; 

Insert a new vertex into c, fully connected to the vertices composing c ; 
Add to C the d + 1 new (d+l)-cliques created by the insertion of the new 
vertex and update r ; 

t «- t+ 1 ; 
end 

Algorithm 2: P-RAN constructing algorithm 

Initially, a network containing d + 2 vertices and d + 2 (d+l)-cliques is created. At 
each time step, a new vertex is added into a (d+l)-clique selected at random. The new 
vertex is connected to each vertex of the selected clique, creating d+1 new (d+l)-cliques. 
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Thus, comparatively to RAN for which new vertices are inserted into (d+l)-cliques that 
contain no vertex inside, for P-RAN, any (d+l)-clique can be selected to insert a new 
vertex, what ever the number of inside vertices is. Figure 13.11 shows the first four steps 
of construction of a P-RAN. A parallel embranchment is created at the third step since 
a clique containing already a vertex is selected for the insertion of a second inner vertex. 



4. Parallel degree distribution and parallel coefficient 

The parallel degree is a characteristic that applies to (d+l)-cliques. We show hereinafter 
that the discrete parallel distribution for P-RANs follows asymptotically a geometrical 
law. 

Definition We define the parallel degree of a (d+l)-clique as the number of vertices 
"inside" the (d+l)-clique, e.g. the number of vertices that are connected to every 
vertices of the (d+l)-clique but are not in the set of vertices that compose the (d+1)- 
clique. 



4-1. Estimating the parallel degree distribution 

Lemma 4.1 Let t be the iteration step of the construction of the growing P-RAN 
algorithm, and let m be an integer. For large t the parallel degree distribution of 
a d dimensional P-RAN asymptotically follows the geometric distribution Pc(m) = 

d+l 



(d+2)(™+ 1 ) ' 

Proof At time t = 0, the networks is composed with d + 2 vertices forming d + 2 
(d+l)-cliques. Each time a new vertex is inserted into the network, the number of 
(d+l)-cliques increases by d + 1. If Nc t is the number of (d+l)-cliques at time t, we 
have Nc t = d + 2 + t.(d + 1). 

Furthermore, each time a (d+l)-clique Cj is selected for the insertion of a new 
vertex, its parallel degree mj increases by 1. Thus, if Nct(m) is the number of (d+1)- 
cliques having a parallel degree equal to m at time t we get the following growth rate 
for Nc t (m) 

at f \ at f \ , Nct-^m-l) Nct_x{m) 

Nct(m) = Nct-i(m) + — ; — Hi) 

v ; u ; d + 2 + (d + l)(t-l) d + 2 + (d+ l)(t- l) v ; 

Let Pc t (m) be the probability to select a (d+l)-clique with parallel degree m at time 
t. Pctijn) can be approximated by the ratio ■ Thus Nc t (m) = Nc t .Pc t (m) = 

(d + 2 + t.(d + l)).Pc t (m) and we get from Eq{T] 



(d + 2 + t.(d + l)).Pct 



m) 



(d + 2 + (t-l).(d+l)).Pct_i(m) 
+Pc t -i{m - 1) - Pc t -x(m) 



Thus 



Pct(m) 



t(d + l) 



d + 2 + t(d + l 



-.Pc t _i(m) + 



Pc t -\(m - 1) 
d + 2 + t{d+ V 



(2) 



(3) 
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As Pc t (m) is bounded for all m and t, from EqJ3]we get that Pc t (m) is a Cauchy 
sequence, which shows that lim^ +00 Pc t (m) = Pc(m) exists and that for large t, 
Pct(m) ~ Pct-i(m) ~ Pc(m). Rewriting the previous equation for large t we get 

Pc{m - 1) Pc(0) 



Pc{m) 



al + 2 



(4) 



(d + 2) 

It is easy to show by induction on t that the probability to select at any time t a 
(d+l)-clique having a null parallel degree is Pc(0) = (d+ l)/(d+ 2). Thus for large t 

d 



Pc(m) 



(d + 2] 



1 

(m+i) 



(5) 



This ends the proof and shows that the parallel degree for P-RAN scales as a geometrical 
distribution. | 



Figure 14.11 gives the parallel degree distribution for P-RANs estimated 
experimentally for each dimension from the construction of 10 networks utterances 
containing 100000 vertices each. The figure gives also the absolute error and 
its corresponding standard deviation measured comparatively to the theoretical 
expectation, showing a good match between simulation and the theoretical model. 



0,8 



■0,001 




Figure 4.1. Parallel Degree distribution estimated from 10 2-dimensional P-RANs 
containing 100000 vertices each. Error and standard deviation to theory are given on 
the right vertical axis. 



4-2. Average parallel degree and parallel coefficient 

Definition The average parallel degree M of a P-RAN is defined as the mathematical 
expectation of the parallel degree, i.e. 

00 °° d + 1 1 

M = E{Pc{m)) = m.Pc(m) = £ m. m+1) = — (6) 

m=l m=l 
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Thus, M measures the average number of vertices inside (d+l)-cliques of a P-RAN. 

Definition We define the parallel coefficient of a d-dimensional P-RAN p as M — Pc(l), 
i.e. 



P 



m=2 



d + 1 



(d + 2)( m + 1 ) 



(7) 



Lemma 4.2 For d-dimensional P-RAN the average parallel degree is M = l/(d + 1), 
and the parallel coefficient is p = (^wf^p • 

Proof According to Eqj6] The parallel degree distribution follows a geometrical law 
whose expectation is M = l/(d + 1) and variance is (d + 2)/(d + l) 2 . Thus p = 
M - Pc(l) = ^ - and the result follows. | 

For d = 2 we get M = 1/3 for P-RAN, which is also the case for RAN, and p = 7/48 
for P-RAN while p = for RAN. 

Figure 14.21 shows the parallel coefficients for P-RANs estimated experimentally 
for each dimension from the construction of 10 networks utterances containing 100000 
vertices each. The figure gives also the absolute error and its corresponding standard 
deviation measured comparatively to the theoretical expectation. 



■0,001 




Figure 4.2. Parallel coefficient of a P-RAN as a function of the dimension. Error and 
standard deviation to theory are given on the right vertical axis. 



5. Estimating the degree distribution 

The degree of a vertex in a network is the number of connections it shares with other 
vertices and the degree distribution is the probability distribution of these degrees over 
the whole network. 
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5.1. Determining the degree distribution 

Lemma 5.1 The degree distribution of a d-dimensional P-RAN is given by the following 
recursion 

' P(k) ~ iEfeSi • P(k - 1) fork>d+l 
P(k) ~ \ fork = d+l K J 



Proof We note that, once a new vertex is added into the P-RAN network, the number 
of (d+l)-cliques available for the insertion of a new vertex is increased by d + 1. After 
t iterations, the number of (d+l)-cliques available for the insertion of a new vertex is 
d + 2 + t(d + l). 

Thus, given a vertex u i; when its degree increases by 1 the number of (d+l)-cliques 
that contain vertex Vi increases by d. So the number of (d+l)-cliques available for 
selection containing vertex Vi with degree ki is (ki — (d + l)).d + d + 1 = d.ki — d 2 + 1, 
since at t — tj, the creation time of vertex Vi there is d + 1 (d+l)-cliques that contain 
vertex Vi. 

Let N t be the total number of vertices into the P-RAN at step t (N t = d + 2 + t) 
and let N t (k) be the number of vertices having a degree k at time t. We can write the 
following difference equation 



N t (k) = N t ^(k) +^S^)N t ^k-l) 

d.k-d 2 +l at (h.\ 
d+2+(t-l).(d+l) iV *-H' L J 



(9) 



Let Pt{k) be the probability to select a vertex with degree k at time t. Pt{k) can 

N t (k) 
d+2+t ' 



be approximated by the ratio j^ff _ Hence Nc t {k) = (d + 2 + t).P t {k) and we get from 



EqJH] 



P t {k).{d + 2 + t)= P t _i(fc).(d + 2 + (t - 1)) 



+lS^mk- P ^( k ~ + 2 + (t - 1)) (10) 
- , +2 f(7-W) ^-i(A:).(d + 2 + (t-l)) 

As -P*(A;) is bounded for all /c and t from Eq{10] we get that Pt{k) is a Cauchy 
sequence, showing that \im t ^ +00 Pt(k) = P{k) exists, and that for large t, Pt{k) ~ 
Pt-i(k) ~ P{k)- Rewriting the previous equation for large t we get 



P(k). (l + d - k ^l +1 ) - (Lk -f + - d+1 .P(k - 1) (11) 
and finally 

f ' t »~ ":^!^ f ' t -" (12) 

This recursive equation is defined for k > d+1. We show next that P(d+ 1) = 1/2 
for all dimensions. 
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• Let Nd+i t t be the expected number of vertices into the network having a degree 
equal to d + 1 at time t, 

• let n t be the expected total number of (d+l)-cliques having a parallel degree equal 
to 0, 

• let n' t be the expected total number of (d+l)-cliques having a parallel degree equal 
to for which all vertices have a degree k > (d + 1) at time t, 

• let n" be the expected total number of (d+l)-cliques having a parallel degree equal 
to for which all vertices have a degree k > d + 1 except one vertex that has a 
degree k — d + 1 at time t. 

For sufficiently large t, every vertex v i in the network has a degree k{ > d + 1, and 
every (d+l)-clique has either all its vertices with a degree k > d + 1 or only one 
vertex with a degree /c = 3. Thus, when we insert a new vertex in a (d+l)-clique Cj, 
only three cases arise for the (d+l)-clique selected for the insertion: 

(i) If the clique has a parallel degree m > then A^ +1 (t) is increased by one, n' t is 
unchanged and n" is increased by d + 1 

(ii) If the clique has a parallel degree m = and all its d + 1 vertices has a degree 
k > d + 1, in this case the N d+ i(t) is increased by one, n£ is decreased by one and 
n" is increased by d + 1 

(iii) If the clique has a parallel degree m = and all its d + 1 vertices have a degree 

> d + 1 except one with a degree equal to d + 1, Nd + i(t) is unchanged, n[ is 
increased by d and is unchanged. 

In Section H] we have shown that the probability to select randomly a (d+l)-clique 
with a parallel degree m = is P(m = 0) = (c?+ l)/(d+ 2) and n t ~ t. ^j . Previous 
statements lead to the following equations 



Assuming that \ivn t _> +00 n' t /n t exists (this is obviously the case since P(k — d+ 1) 
exists), n' t ~ a.t where a is a constant. Replacing n' t in Eq{13]we get 



P(d+1) 



1 d+1 n' t 



(13) 



d + 2 d + 2'n t 




(14) 




(15) 



leading to a = f. Thus, 
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Finally, P(d + 1) = P(k = d + 1) ~ ^ + J = g + ^1) = 1/2. Note that 
P(g? + 1) is independent from the dimension d. 

This completes the recursive equation that gives the degree spectrum distribution. 

I 

To our knowledge, there is no simple analytical expression for P(k) in any 
dimension. Nevertheless, for d — 1, we get 

This result in dimension one has already been reported in [5]. 
5.#. Degree exponent 

For scale free networks, the degree distribution follows asymptotically a power law whose 
exponent is called the degree exponent. In the following, we show that P-RANs are scale 
free networks and derive their degree exponents. 



Lemma 5.2 The degree exponent of a d- dimensional P-RAN is 7 



2.d+l 



Proof To show that the degree distribution follows a power law, we evaluate the 
asymptotic value of the following ratio 

P (b] _ l °9(P(k)) - log{P{k - 1)) _ log(P(k)/P(k - 1)) 
{ ) log(k) - log(k - 1) log(k/(k-l)) K ' 

Thus 



R(k) 



ay d.k-d 2 +d+2' " y d.k-d 2 +d+2 ' 

log(k/(k-l)) ~ log(k/(k-l)) 

-d 2 



A+- 



=ggjj- ) 



(19) 



1+ 



-log{l-\) 



and for large k 

R{k) ~k.( 



-d 2 -d+l -d 2 +d+2 



d.k d.k 

2.d+l 



(20) 



d 

This shows that for large k P{k) ~ k~~< with 7 = 2,d , +1 . 



For d = 2, we theoretically get 7 = 5/2. 

We evaluate the empirical degree exponent using the mean of the maximum 
likelihood estimate computed according to the following formula proposed in [3] : 
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Xmin 

Figure 5.1. Degree exponent estimation for a 2-dimensional P-RAN containing 
200000 vertices according to Eq|2i1 

where fcj, i = l,2,...,n are the observed values of k such that ki > k min . Figure I5TT1 
gives the estimated degree exponent according to k min for networks containing 500000 
nodes. Results are well correlated with theory for k min £ [30; 300]. When x m i n is lower 
than 30, a bias is introduced in the power law by low degree vertices while when k min is 
higher than 300, the set of vertices having a high degree becomes too small to give an 
accurate estimate. 

6. Clustering coefficients 

The clustering coefficient Ci that characterizes vertex Vi is the proportion of links 
between the vertices within its neighborhood (vi excluded) divided by the number of 
links that could possibly exist between them. For undirected graph, considering two 
vertices Vi and Vj, the edges Vi — > Vj and Vj — > are considered identical. Therefore, 
if a vertex v i has ki neighbors, k ^ kl ~ l > edges could exist among the vertices within its 
neighborhood. The clustering coefficient for the whole network is the average of the 
clustering coefficients Ci over the set of vertices composing the network, i.e. this is the 
expectation of the clustering coefficient distribution. 

When a vertex is inserted into the network, it is connected to all the vertices of 
a selected (d+l)-clique. It follows that every vertex having a degree ki — d + 1 has a 
clustering coefficient of one. Furthermore, when a vertex Vi having a degree ki belongs 
to a (d+l)-clique in which a new vertex is inserted, its degree increases by one and 
the new inserted neighbor connects to d vertices among the ki vertices that compose 
its neighborhood previously to the insertion. This leads to the following clustering 
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coefficient for a vertex having a degree k 



C(k) 



d.(d+l) 
2 



+ d.(k-d-l) d.(2k-d-l) 



fc.(fe-i) 



k.{k-l) 



(22) 



This local clustering coefficient is exactly the same as the one obtained for vertices 
in RAN [13] . Eqf22l shows that the local clustering coefficient scales as C(k) ~ k^ 1 

We average these coefficients using the discrete degree distribution (EqJSJ) as follows 



C 



d.(2ki -d-1] 
k{.{ki 1) 



■P(ki 



(23) 



ki=d+l 

For d — 2, we get C = 0.813. Figure 16.11 shows that the clustering coefficient 
increases from 0.813 for d = 2 to 1 as d tends towards infinity. Comparatively, HDRANs 
have a significantly lower clustering coefficient at low dimension, e.g. for d — 2, a RAN 
has a clustering coefficient C = .768. 

Figure 16.11 gives the clustering coefficients for P-RANs estimated experimentally 
for each dimension from the construction of 10 networks utterances containing 100000 
vertices each. The figure gives also the absolute error and its corresponding standard 
deviation measured comparatively to the theoretical expectation, showing a good match 
between simulation and the theoretical model. 



-0,004 




Dimension 

— CC experimental -•- Error with theory -±- Standard deviation 



Figure 6.1. Clustering coefficient of a P-RAN as a function of the dimension. Error 
and standard deviation to theory are given on the right vertical axis. 



7. Average path length 

The average path length (APL) is a characteristic of the network topology that is defined 
as the average number of edges along the shortest paths for all possible pairs of network 
vertices. Following exactly the derivations already presented in [TBI [Ti] for RAN, we 
address below the average path length for P-RAN. 
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First, we suppose that any vertex of the P-RAN network is ordered according to 
its insertion time stamp t that we consider discrete (t e N). It is straightforward to 
establish that for P-RAN the following property holds (as well as for DAN or RAN) 

For any two arbitrary vertices i and j all shortest paths from i to j does not pass 
through a vertex k if k > max(i,j). 

Let d(i,j) denotes the distance between vertices i and j, namely the length of a 
shortest path between vertices i and j. Let c(N) be the sum of all distances between 
all the pairs of vertices into the network with order N, e.g. containing N vertices. 

a(N)= ( 24 ) 

l<i<j<N 

and let L(N) be the average path length of the P-RAN of order N 

Following exactly the approach given in [13J we get the following recursive inequality 
for a(N) 

a(N + 1) < a{N) + N + (26) 

Considering the inequality EqfSn] as an equation we get the same upper bound for 
the variation of c(N) than for RAN 

da(N) Ar 2a(N) 
which leads to 

a(N) < N 2 .log(N) + S (28) 

where S is a constant. As o~(N) is asymptotically upper bounded by ~ N 2 .log(N), 
L(N) is asymptotically upper bounded by log(N), e.g. L(N) increases at most as log(N) 
with N. 

Figure 17.11 compares for dimensions 2, 4 and 6 average path lengths for HDRANs 
and P-RANs and shows that, for a given dimension, average path lengths are shorter 
for P-RANs than for HDRANs. Nevertheless, as the dimension increases, differences 
between path lengths vanish. This result was expected since P-RANs have a higher 
clustering coefficient than RANs. 



8. Conclusion 



From previous works on Apollonian Networks, mainly RAN and HDRAN networks, 
we have introduced what we call Parallel Deterministic or Parallel Random Apollonian 
Networks. These topologies, for which (d+l)-cliques may accept in their inside more 
than one vertex fully connected to the vertices composing the clique, are still small world 
and scale free. This paper reports the main statistical properties of P-RANs. For such 
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Figure 7.1. Average path length in RANs and P-RANs 

networks, the degree exponent is in between 2 (2 being attained at the limit when the 
dimension tends towards infinity) and 2.5 (when the dimension of the network is 2) or 
3 if we accept the limit case of Apollonian networks in dimension one. We have shown 
analytically that, comparatively to RAN or HDRAN, P-RAN networks are characterized 
with higher clustering coefficients and shorter average path lengths. P-RAN are also 
characterized by their parallel degree distribution and parallel coefficient that quantify 
the number of vertices inside the (d+l)-cliques that compose P-RAN networks. The 
simulations results provided through out the paper are in very good conformance with 
the analytical expectations. 
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